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I. Introduction 


In order to understand the physical and chemical processes which produce the observed 
spatial morphology of the cometary coma, it is necessary to analyze observational data with 
physically meaningful models. With this in mind, the goals of this project have been to undertake 
a coupled program of theoretical modeling and complementary observational data analysis 
regarding the spatial distributions of neutral gases in the coma. More specifically, the particular 
topics of interest are: (1) the theoretical modeling of the non-equilibrium dynamics of the inner 
coma with emphasis on the region of the coma from the transition from collisional fluid flow out to 
the free-flow region and on observable conditions in the coma (i.e. density, outflow speed, and 
temperature), and (2) the model analysis of an important set of long-slit CCD spectra of comets. 

A range of approaches for modeling the cometary coma have been employed by many 
investigators. The Monte Carlo and the coupled hydrodynamic/Monte Carlo approaches to 
modeling the cometary coma recognize that the region of the coma typically observed is that where 
collisions do occur between molecules, but where not enough collisions occur to be able to 
describe the coma as a fluid. Furthermore, outside this transition region, the Monte Carlo model 
have already been demonstrated to be a powerful tool for representing the complexities of free 
molecular flow, such as time-dependence, radiation pressure, the orbital mechanics of individual 
particles, etc. The theoretical modeling involves the continued development and use of Monte 
Carlo models and coupled hydrodynamic/Monte Carlo models for the outflowing cometary coma. 

The observational data analysis portion of this project is an integral and important part of the 
overall picture. Developing sophisticated models is not an end unto itself. The models are only 
important if they can explain and be tested by observational data. Conversely, models are used for 
the interpretation and analysis of observational data by necessity. The side-by-side development of 
models along with the observation and analysis of data is an important and integral part of this 
project. The understanding of data from Halley is a valuable endeavor for the preparation of 
NASA’s CRAF (Comet Rendezvous Asteroid Flyby ) mission. The scientific community has in 
hand valuable observational and in situ data regarding one comet (Halley). It is important to use 
Halley as the benchmark by which other remotely observed comet data can be understood. The 
self-consistent analysis of data with appropriate models is therefore of the utmost importance. 

The data analysis work for this project includes the analysis of the spatial profiles of [01], 
NH2, CN and C2- It is being carried out in collaboration with Dr. Uwe Fink of the University of 
Arizona who with his co-workers observed comet Halley extensively during both the pre and post 
perihelion periods. 
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II. Progress during this Year 


Significant progress was made this year in both the theoretical and the observational data 
analysis areas. Each area will be discussed separately. 


II.A. Theoretical Modeling 


In a collisional Monte Carlo model for the cometary coma, the trajectories of molecules are 
calculated explicitly, following the the complete time history of the particle from its production at 
the nucleus through various stages of photodissociation, entrainment within the outflowing coma, 
and/or inter-molecular collisions. Between two collisions, the path length of a particle is calculated 
given a set of random numbers and a mathematical description of the conditions regarding the state 
of the background gas through which it moves. These conditions are the temperature and bulk 
motion of the gas. Two special cases of the collision path length were presented in the paper by 
Combi and Smyth (1988a). One was for the important limit which is correct for the fast moving H 
atoms and yields a collision path integral which is integrable and therefore fast to compute in a 
model. It essentially assumes that the coma is stationary and has negligible temperature. The other 
case was the more general one which accounted for the relative motion of the particle with respect 
to the outflowing coma, but where the coma temperature was still neglected. It should be noted 
here that for the purpose of computing gas density (or similarly the resulting column density) 
profiles, the simpler formulations of Combi and Smyth (1988a) are quite adequate. 

On the other hand if one is interested in information about the velocity distribution function, 
then neither of these formulations treats nearly thermalized heavy molecules very well. Therefore, 
we have performed the derivation of the more general case of the collision rate at any location in a 
radially outflowing coma with a given Maxwell-Boltzmann temperature. It has been incorporated it 
into the model. The collision rate for a particle in the coma at a distance r Q from the center of the 
nucleus of the comet is obtained by integrating the relative velocity weighted velocity distribution 
function over the three-space velocity. It is given by 


v=-j^J~[erf(U^){l + 
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c 
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exp(-U_ a) 

1 Vi aU 1 


where 


2 


U = (u 2 -2uc 0 cos9 + c Q 2 ) U2 
P = c 0 lu 
a = m / 2kT 


Co = coma outflow speed 
u = particle speed 

0 = particle direction angle relative to the outward radial 
= collision cross section between the particle and the coma gas 
Q = coma gas production rate. 


In the Monte Carlo model a random number on the interval from 0 to 1 is set equal to the 
collision probability along the path and the collision path length corresponding to that random 
number is then chosen. Since the outflow speed and temperature vary significantly only for large 
changes in radial distance, we can assume (as in Combi and Smyth 1988a) that over any one 
collision path they are both constant. This calculation results in the following integral equation: 


-ln(l-R) = ~T^~ j 0 
4nc 0 r o y J o 


1 expf-v,. 2 ) 

v x .lerf(v x .){l + — ; + ** X — ] 

Q° 2 V dX > 


1 ~f~ 2X' cos Q -h X' 


where 


.. f o X + cos 0 . n 2 i 11 2 

v =y[ 1-2 p 2 ' I T T + P 1 

A (1 + 2 X cos 0+ X ) 


7 = u^y~a 

A' = collision path length normalized by tq. 


Unlike the special case for fast H atoms this integral cannot be solved analytically and 
inverted. However, we have developed an approximation scheme which yields a single-iteration 
solution, based on the fact that in realistic cases the numerator within the integral is only a weakly 
varying function along the path length. The approximation is a two-step predictor corrector type 
solution. The first step is to make an initial guess whereby we assume that the numerator in the 
integral is a constant. This accounts for the temperature and relative velocity of the particle near the 
start of the collision path. With this term as a constant the remaining integral can be done 
analytically as in Combi and Smyth (1988a). This path length serves as a first guess or predictor 
of the true path length. As a second step we then replace the numerator of the real integral with its 
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value at the half-way point along the predicted collision path and proceed to recalculate the a 
corrected collision path. 

This was tested by calculating the collision probability along the path from some point of 
origin out to infinity both with the new method and by explicit numerical integration. We 
concentrated on the expected troublesome cases where the method might be expected to do the 
worst. In the inner part of the coma where the collision paths are all short because the densities are 
high the method gives almost an exact result. Fortunately in realistic models there are very few 
particles which have velocity vectors in the regime where the approximation is not as good. 

An example of a pathological case would be a radical which is located at 2-5 collision path 
radii from the nucleus and which has a velocity magnitude equal to the outflow speed but is 
directed to within an angle of 10 to 30 degrees of the negative radial direction. Note that it is very 
improbable to generate such radicals as dissociation products of outflowing parents. However, 
even in this case it is not until the collision probability exceeds 95% that the approximate collision 
path length begins to diverge only somewhat from an actual collision path integral calculation. The 
conclusion is thus that only a very small fraction of collisions for only a small fraction of molecules 
will have some small error in the calculation of collision paths. 

Once a collision occurs, we choose a randomized target particle from the background gas, 
which is a moving Maxwellian. Since it is easier to choose the background gas particle from an 
isotropic (in angle) moving Maxwellian rather than from the correct relative velocity weighted 
distribution, we have employed the principle of the 'prejudiced' source - or the skewed probability 
density distribution of Combi and Smyth (1988a) - where we choose from one distribution and 
then weight the particle according to the relative magnitude of the real to the assumed probability 
density distribution. This prevents us from needing to generate look-up tables for complicated 
multidimensional functions. This also minimizes both the computation times for model runs and 
development time in changing the model. These new features have been incorporated into the 
model. These enhancements improve our modeling of the typical cometary radicals and represent 
an important step toward studying the fundamental kinetic nature of much of the coma. 

We are now planning to apply this model to realistic cases of NH2, C2 and CN from the 
analysis portion of the project during the continuation grant of this project. 

II.B. Observed Spatial Profile Analysis 

Six months into this current project year our goal was the completion of the work and the 
preparation of a paper which deals with analysis of a number of spatial profiles of [01] and NH2 
determined from the long-slit CCD spectra of Comet Halley taken by Drs. Uwe Fink, M. DiSanti 
and A. Schultz at the University of Arizona. A summary of the data appears in Table 1. 
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Table 1. Summary of Long Slit Spectroscopy of Comet P/Halley 


Date 

1985-86 

r 

(AU) 

A 

(AU) 

Number of Spectra 

Aug 23 

2.84 

3.24 

3 

Aug 24 

2.83 

3.21 

1 

Sep 17 

2.53 

2.49 

3 

Sep 24 

2.43 

2.27 

3 

Oct 20 

2.09 

1.43 

5 

Oct 21 

2.07 

1.40 

2 

Nov 15 

1.72 

0.74 

9 

Dec 8 

1.38 

0.70 

10 

Dec 9 

1.36 

0.72 

5 

Jan 10 

0.87 

1.32 

13 

Jan 11 

0.86 

1.34 

11 

Jan 12 

0.84 

1.36 

9 

Marl 

0.72 

1.27 

15 

Mar 2 

0.74 

1.25 

11 

Mar 3 

0.75 

1.22 

6 

Mar 4 

0.76 

1.20 

6 

Mar 5 

0.77 

1.18 

9 

Apr 14 

1.38 

0.43 

22 

Apr 15 

1.39 

0.44 

19 

May 9 

1.75 

1.05 

14 

Jun 5 

2.13 

1.93 

6 

Jun 6 

2.14 

1.96 

5 


r = heliocentric distance in AU 
A = geocentric distance in AU 
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Spatial profiles of cometary C 2 and CN were to be produced by Dr. Fink and co-workers 
for a subsequent paper during the following year (during the next grant period). When we began 
looking at the spatial profiles from the spectra, it became evident that sky background subtraction 
procedure adopted by Fink and co-workers was effecting the shape of the spatial profiles far from 
the nucleus. It is the level at distant locations along the profile which is the most important from 
the standpoint of inferring lifetimes and scale lengths and constraining model parameters. 

Because of these complications Dr. Fink has revised his procedure for performing the sky 
level assessment and has re-reduced the spatial profiles for all four species from most of the CCD 
spectra. The result of this has been to push back the publication date of the first paper, which will 
deal with a simple analysis of the spatial profiles of all four species, until early in the next project 
year. Whereas, the publication of the second paper, which will deal with a more substantial model 
analysis concentrating on C 2 and CN will still be finished on the original schedule by the end of 
the first project year of the follow-up grant. 

By the end of this current project year we have received all of the revised spatial profiles for 
all four species and have performed the first-step Haser model analysis. There are some 
restrictions regarding the fitted scale lengths and the fitting procedure, which will be discussed here 
species by species. An Appendix to this report contains all of the observed profiles with best fit 
Haser models. 

II.B.a. [01] 

The forbidden oxygen emission is now believed to originate from oxygen atoms produced in 
the metastable ^D state upon photodissociation of their parent molecules. The major sources of 
O(lD) atoms are now believed to be H 2 O in the inner coma (Biermann and Trefftz 1964) and OH 
in the outer coma (van Dishoeck and Dalgamo 1984). Roughly 100 seconds after the atoms are 
produced they emit either the 6300 A or the 6364 A photon (in the ratio of 3/1), unless they are 
produced in the very inner coma (<1000 km) where the state can be collisionally quenched 
before a photon emission (Festou and Feldman 1984). The result of this is that the 0(1 D) 
emission is proportional to the H 2 O column in the inner coma (<10^ km) and to the OH column in 
the outer coma (>10^ km). Magee-Sauer et al. (1989) have published simple model profiles for 
the combined source given alternative sets of branching ratios. 

The end result is that we expect the spatial distribution of 0(1 D) emission to decrease as the 
inverse of the distance to the nucleus and have a scale length for decay somewhere intermediate 
between the scale lengths for decay of H 2 O and OH. Table 2 gives the results for the O(lD) 
profiles. In fact all except one of the spatial profiles is well reproduced by a point source model 
with a decay scale length, except for the case of December profile where a still very small source 
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Table 2. Haser Scale Lengths for [01] Profiles 



observed 

m 2 

m \.5 

Date 

P D 

P D 

P D 


Oct 

0 

49 

0 

1 1 

0 

16 

Dec 

2.1 

330 

1.1 

175 

1.3 

205 

Jan 10 

0 

40 

0 

54 

0 

50 

Jan 11 

0 

54 

0 

77 

0 

70 

Jan 11C 

0 

330 

0 

446 

0 

428 

Jan 12 

0 

48 

0 

65 

0 

62 

Marl 

0 


0 


0 


Mar2 

0 

33 

0 

60 

0 

52 

May 

0 

- 

0 

- 

0 

“ 


P = Parent scale length in 1C)3 km. 

D = Daughter scale length in 1(P km. 
rjp^ = Reduced to 1 AU by an rp^ law 
rf^l-5 = Reduced to 1 AU by an ipp^^ law 
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scale length is required for a model fit. This is likely, however, to be a fictitious scale length 
whereby the flattening is due simply to the irregular activity of the gas production. 

Unfortunately, except for the December and the January 1 1 composite profile (created from a 
regular exposure and one with a long integration time where the nucleus was masked), it is 
possible that we are not really seeing the true effective decay scale length but simply the fact that 
the profile becomes too noisy. before a significant departure from 1/r can be seen. The fact that we 
do see the 1 /r distribution near the nucleus implies that the reduction procedure, the subtraction of 
of the contribution to the 6300 A emission from NH 2 and the evaluation of the sky level are all 

reasonably well done. 

II.B.b. NH2 

For NH2 there just are not many published spatial profiles and determined scale lengths 
available. It has been suggested in recent, but as yet unpublished results, that the NH 2 spatial 
distribution is consistent with primary production from the photodissociation of NH 3 (Wyckoff et 
al. 1990) and with decay primarily through photodissociation producing the visible NH (Schleicher 
and Millis 1989). Table 3 gives the results for the NH 2 Haser model scale lengths. When reduced 
to 1 AU by an rn* ^ law. This variation is expected from the combination of photodissociation of 
a parent molecule, whose lifetime varies as rH“> and the approximate rH ' 05 variation in the coma 
outflow speed from photochemical heating (Combi 1989). Since the NH 2 radicals are produced 
generally within 10,000 km from the nucleus they are also expect (at least for the heliocentric 
distances covered) to be thermalized and follow the same outflow speed as the water coma. The 
straightforward power law fit in rn in fact agrees with this interpretation. 

Given a variation of rn^ we find values for the lifetimes of the NH 2 parent and NH 2 , 
respectively to be 5500 km and 40,000 km. Given an outflow speed, at 1000-10,000 km from the 
nucleus of about 0.8 km/s for NH 3 this implies a lifetime for NH 3 of 6900 seconds which agrees 
quite well with the expected value of 7700 seconds (Allen, 1989 private communication). A 
lifetime for NH 2 with an average velocity closer to 1.0 km/s would yield a lifetime of about 40,000 
seconds at 1 AU. Note that for immediate (or nearly immediate) thermalization, the radial motion 
assumed by the Haser model is quite realistic and and appropriate to first order. OH, on the other 
hand, is produced much farther from the nucleus and vectorial ejection of the the radicals must be 
taken into account (Festou 1981). We will test the thermalization of NH 2 , in the context of these 
data, using the new heavy molecule collision algorithm discussed in the previous section of this 
report. 
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Table 3. Haser Scale Lengths for NH2 Profiles 



observed 

r H 2 

rH 1 ' 5 

Date 

P D 

P D 

P D 


Oct 

12.2 

365 

2.8 

84 

4.1 

122 

Dec 

9.4 

109 

5.0 

58 

5.9 

68 

Jan 10 

4.0 

46.5 

5.3 

61 

4.9 

57 

Jan 11 

2.8 

59.4 

3.4 

80 

3.5 

74 

JanllC 

2.6 

72.3 

3.5 

98 

3.3 

91 

Jan 12 

2.6 

63.2 

3.7 

90 

3.4 

82 

Marl 

4.4 

11.9 

8.5 

23 

7.2 

19 

Mar2 

4.6 

13.7 

8.4 

25 

7.2 

22 

May 

21.4 

52.9 

7.0 

17 

9.2 

23 


P = Parent scale length in 10^ km. 

D = Daughter scale length in 10^ km. 
rn 2 = Reduced to 1 AU by an rn 2 law 
rH^ = Reduced to 1 AU by an ^ law 
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II.B.c. CN 

The last summary of pre-Halley spatial profiles for CN was done by Combi and Delsemme 
(1986). In that paper we found that from data for a number of comets that the parent Haser scale 
length was about 2 x lO^ km and seemed to vary as rn^-^ which is (like NH2 above) consistent 
with production from the photodissociation of a parent molecule. The daughter scale length 
seemed to be of order 3.2 x 10^ km at 1 AU but the heliocentric distance dependence was always 
too irregular to determine. For the new spatial profiles from comet Halley we find that some of the 
profiles do not extend to large enough distances from the nucleus to enable the determination of 
both scale lengths. For these we assume a value consistent with the value at 1 AU which scales 
with a square law dependence. Since the daughter scale length is very long compared with the 
parent value, the parent scale length determined this way is only weakly dependent on the assumed 
daughter value anyway. For some we can and have determined both scale lengths. Table 4 gives 
the list of best fit Haser models for the CN profiles. 

The values determined so far are reasonably consistent with our pre-Halley summary, 
although we are now in the process of assembling the results from a number of investigations in 
which CN scale lengths were determined for Halley. These results will be compiled with our new 
data results in the second paper to be finished by the end of the first year of the continuation project 
to this one (Combi, Fink and Di Santi 1991, in preparation). 

II.B.d. C2 

Perhaps the most surprising results come from the Haser model analysis of the C2 profiles. 
In the pre-Halley summary (Combi and Delsemme 1986), we found that both the parent 1.6 x 10^ 
km) and daughter (1.1 x 10^ km) scale lengths varied as the square of the heliocentric distance. In 
this new data we have profiles which extend both very close to the nucleus and fairly far away, as 
compared generally with the data we had in our 1986 summary. In the old data the profiles never 
did go as close to the nucleus. 

For all of the new C2 profiles, we require a Haser model where the parent and daughter scale 
lengths are equal. Furthermore, the value varies as the square of the heliocentric distance as 
opposed to the slower variation implied by the CN, and NH2- Note that for the old summary 
Combi and Delsemme (1986) found a ratio of about 7 between the daughter and the parent. The 
reason an equal scale length model is required is because all of the profiles are very flat near the 
nucleus and turn over very rapidly. This behavior has been noted by other investigators (Cochran 
1985; Wyckoff et al. 1988) going as far back as Delsemme and Miller (1971). 

Various suggestions for the unusual appearance of the C2 profiles including a 2-step 
photodissociation (grandparent-to-parent-C2 ) and a grain source have been made (see Wyckoff et 


Table 4. Haser Scale Lengths for CN Profiles 



observed 

r H 

2 

rH 

1.5 

Date 

P 

D 

P 

D 

P 

D 

Oct 

- 

- 

- 

- 

- 

- 

Dec 

24 

550 

12.4 

293 

14.9 

343 

Jan 10 

22 

413 

29 

546 

27.1 

509 

Jan 11 

25 

376 

34 

508 

31 

471 

Jan 11C 

31 

(400) 

23 

(540) 

39 

(501) 

Jan 12 

21 

(375) 

14.8 

(531) 

27 

(487) 


Marl 

15 

730 

29 

1408 

25 

1195 

Mar2 

45 

45 

82 

82 

71 

71 

May 

44 

(980) 

14 

(320) 

19 

138 


P = Parent scale length in 10^ km. 

D = Daughter scale length in 1(P km. 

Tfp. - Reduced to 1 AU by an rpj^ law 
rf|L5 = Reduced to 1 AU by an rn^5 law 

Values in parentheses are daughter scale lengths assumed either from the value on an adjacent day 
or from the 3.2 x 10^ km extrapolated to the actual heliocentric distance. 



al. 1988). Continued work on the unusual appearance and behavior of the C 2 profiles and their 
heliocentric distance dependence will be one of the principal subjects of the second paper to be 
written later this year (Combi. Fink and Di Santi, 1991, in preparation). Table 5 shows the list of 
the best fit Haser scale lengths. 


Table 5. Haser Scale Lengths for C 2 Profiles 


Date 

Observed 

7 

ITT 

Oct 

255 

59 

Dec 

107 

57 

JanlO 

40 

53 

Jan 11 

41 

55 

Janl 1C 

44 

59 

Jan 12 

39 

55 


Marl 

26 

50 

Mar2 

28 

51 

May 

284 

93 


= Reduced to 1 AU by an rf-j- law 

For all C 2 profiles the parent and daughter scale lengths are equal and in units of 
10 3 km. 
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APPENDIX 


A Comparison of Haser Models vs. Observed Profiles 

The following is a set of Figures for all of the determined spatial profiles of O(^D), NH 2 , 

CN and C2- Each figure is a plot of Log Brightness (Arbitrary Units) vs. Log Distance from the 
Nucleus in (km). For some of the figures the fitted scale lengths are printed at the bottom. Below 
is a Table containing the coded name for each of the spatial profiles giving the date and species. 
The coded name contains the name of species also. In some cases the spatial profile was created 
from the sums of spectra taken on consecutive days (e.g. Oct 20-21, Dec 8-9). There is a set of 
profiles taken on Janqary 1 1 (denoted as Jan 1 1C) which was created as a composite of a "normal" 
spectrum plus a "masked" spectrum of long integration time which yields a better signal-to-noise 
ratio than the normal spectra alone. The orbital geometry parameters are given in Table 1 in the 
main body of this report. 


Table A.l. Identification of Spatial Profile Name Codes 


Date 

1985-1986 

O(lD) 

NH2(0,8,0) 

C2(l-0) 

CN(l-O) 

(red) 

Oct 20-21 

i85oicsu 

i85n08sb 

i85c2sum 

i85cnlsu. 

Dec 8-9 

kl90oicd 

k85n08su 

kl90c2cd 

kl90cncd 

Jan 10 

al74oib 

al74nh8 

al74c2 

al74cnl 

Jan 11 

janoiasu 

jann08su 

janc2sum 

jancnlsu 

JanllC 

jlloiasu 

j 1 1 n08cd 

jl lc2cbd 

jl lcnlcd 

Jan 12 

jl2oiasu 

jl2n08su 

jl2c2sum 

jl2cnlsu 

Marl 

mOloibsu 

m01nh8su 

m01c2sum 

mOlcnlsu 

Mar 2 

m02oibsu 

m02nh8su 

m02c2sum 

m02cnlsu 

May 9 

e86oisum 

e86n08su 

e86c2sum 

e86cnlsu 


1 6 
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